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ABSTRACT 


An  analytical  study  is  made  of  the  flow  past  hydro¬ 
foils  with  large  camber  at  arbitrary  angles  of  attack. 
The  class  of  hydrofoils  considered  is  restricted  to 
those  whose  hodograph  plane  has  a  shape  very  near  that 
bounded  by  two  circular  arcs.  Under  this  assumption  a 
general  method  is  derived  to  calculate  the  parameters  of 
the  flow.  In  addition,  a  short-cut  procedure  is  intro¬ 
duced  for  estimating  the  lift  and  drag  coefficients.  As 
a  numerical  example  to  illustrate  the  theory  the  flow 
past  a  circular-arc  foil  is  considered.  The  force 
coefficients  are  in  good  agreement  with  the  values 
computed  by  Rosenhead. 


CONTENTS 


Abstraot  . . ii 

1.  Introduction  ........  .  .  .  ...  .  .  .  .  1 

2.  Hydrofoils  with  Zero  Cavitation  Number  ...  3 

3.  Method  of  Analysis  . .  5 

4.  Hydrofoils  with  finite  Cavitation  Number  .  .  12 

5.  A  Rapid  Method  for  Estimating  Lift  and  Drag 

Coefficients  .  14 

Acknowledgments . l6 

Appendix  A.  Perturbation  of  the  flow  Past  a 

flat  Plate . 1 6 

Appendix  B.  List  of  Symbols . 30 

References  . . 3  2 

Tables . . .  34 

Figures  .  ..... 

iii 


1.  INTRODUCTION.  If  a  well-streamlined  body  is  moved 


through  water  at  a  sufficiently  high  speed,  cavitation 
will  ocour  and  become  more  severe  as  the  speed  increases. 
Beoause  of  the  noise  and  structural  damage  whioh  ordinar¬ 
ily  ensue,  elaborate  precautions  are  often  taken  to  avoid 
cavitation.  Experience  shows,  however,  that  at  the 
design  speeds  proposed  for  present-day  hydrofoil  craft, 
cavitation  is,  practically  speaking,  unavoidable.  More¬ 
over,  if  the  zone  of  cavitation  becomes  so  large  that  a 
large  vapor  cavity  is  formed  and  remains  attached  to  the 
body  ( supercavitating  flow),  the  flow  regime  is  completely 
changed  and  damage  is  often  avoidable.  The  point  of  view 
now  taken  is  to  assume  that  cavitation  must  occur  and  to 
design  hydrofoils,  struts,  and  appendages  accordingly. 

In  regard  to  lifting  hydrofoils  and  propellers,  the 
analytical  work  of  Tulin  (1953)  on  supercavitating  foils 
has  been  of  special  significance  because  of  the  compara¬ 
tive  simplicity  with  which  design  parameters  can  be 
estimated,  provided  that  the  camber  and  attack  angle  are 
small.  The  recent  extension  of  Tulin's  theory  by  Chen 
(1962)  to  include  second-order  terms  has  increased  the 
range  of  usefulness  to  intermediate  values  of  oamber  and 
attack  angle  (up  to  15  or  20  degrees,  say),  Chen's  theory, 
which  again  is  remarkably  devoid  of  computational  diffi¬ 
culties,  should  suffice  for  the  great  bulk  of  practical 
cases.  For  certain  special  configurations,  however,  it 
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will  be  necessary  to  have  in  hand  a  theory  which  applies  to 
arbitrary  foil  shapes  having  large  attack  angle,  large 
camber,  or  a  combination  of  both.  Moreover,  it  is  essential 
to  have  available  a  more  inclusive  theory  against  which  to 
check  the  limits  of  application  of  the  approximations  of 
Tulin  and  Chen. 

If  the  angle  of  attack  is  large  but  the  camber  is  small, 
an  excellent  approximate  theory  might  be  developed  by  per¬ 
turbing  about  the  classical  solution  for  a  flat  plate  at 
arbitrary  incidence.  The  structure  of  this  theory  is 
accordingly  presented  in  Appendix  A;  however,  no  numerical 
results  have  as  yet  been  worked  out. 

To  develop  a  more  general  analytical  method,  one  would 
naturally  turn  to  the  non-linear  theory  of  Levi-Civita 
(presented,  for  example,  by  Milne-lhompson,  i960).  The 
method  has  been  applied  by  Brodetsky  (1923),  Rosenhead 
(1928),  and  Wu  (1956)  to  a  limited  variety  of  foil  sections, 
but  the  complexity  of  the  calculations  leaves  something  to 
be  desired.  More  recently,  Wu  (1962)  has  devised  a  non¬ 
linear  theory  which  overcomes  most  of  the  difficulties  (Wu 
and  Wang,  1963).  However,  in  general,  there  appears  to  be 
justification  for  developing  alternative  lines  of  attack  on 
the  general  problem  for  two  reasons.  First,  each  theory 
will  have  an  inherent  advantage  in  calculating  a  certain 
class  of  foils;  for  example,  the  method  of  Wu  should  be 
especially  powerful  in  dealing  with  hydrofoils  with  flaps. 
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Second,  the  various  schemes  all  being  in  some  sense 
approximate,  a  ehanoe  is  afforded  of  maintaining  more  con¬ 
trol  on  the  aoouraoy  by  ohecking  the  result  of  one  analysis 
against  another. 

In  the  present  work,  therefore,  a  theory  is  developed 
to  deal  with  a  oertain  class  of  hydrofoils  having  rather 
large  camber  and  arbitrary  angle  of  attack.  The  method  of 
analysis  is  such  that  the  estimates  of  lift  and  drag  are 
expected  to  be  very  good  if  the  foil  has  a  smoothly  curving 
wetted  camber  line,  and  a  hodograph  plane  nearly  bounded  by 
two  circular  arcs.  The  efficacy  of  the  method  is  illus¬ 
trated  by  comparison  with  some  numerical  results  given  by 
Rosenhead  (1928).  It  almost  goes  without  saying  that  a 
thorough  review  and  comparison  of  the  several  methods  now 
available  would  be  most  welcome. 

2.  HYDROFOILS  WITH  ZERO  CAVITATION  NUMBER.  Let  a  curved 
hydrofoil  abc  be  set  at  an  angle  a  in  an  infinite 
stream  with  uniform  velocity  U,,,  and  a  pressure  p^  at 
infinity  (Figure  la).  For  sufficiently  high  values  of 
the  flow  will  separate  at  both  edges  of  the  foil  leaving  a 
cavity  inside.  If  the  pressure  inside  pfl  is  equal  to  p^, 
then  the  cavitation  number  will  be  equal  to  zero.  The 
cavitation  number  a  is  defined  by  the  formula 


in  which  p  is  the  fluid  density.  In  this  oase  the  free 
streamlines  ol  and  al  will  extend  to  infinity.  Take 
ox  to  be  the  x-axls  and  the  perpendioular  at  o  as  the 
y-axis.  Tor  plane  irrotational  flow  the  ooaplex  potential 
V  can  be  defined  as 

W  -  <P  +  it  (1) 


Here  <p  is  the  Telocity  potential  and  t  is  the  stream 
funotion.  If  point  b  is  taken  as  the  origin  for  the  V- 
plane,  it  will  be  as  shown  in  Figure  lb.  Ihe  complex 
relocity  u>  is  defined  to  be 


u> 


1  dV  u 

TE35  "  K 


(2) 


Here  u  and  t  are  the  Telooity  components  in  the  x  and 
y  directions  respeotiTely .  Ihe  flow  is  assumed  to  be 
steady  and  the  pressure  in  the  caTity  to  be  oonstant; 
consequently,  the  magnitude  of  the  Telooity  on  the  free 
streamlines  al  and  ol  is  oonstant  and  equal  to  U,*. 

Thus  these  streamlines  are  represented  by  a  portion  of  the 
oiroumferenoe  of  a  unit  oirole  in  the  w-plane  (Figure  lo). 
The  Telocity  components  at  points  o  and  a  are  giTen  by 

2  “I/2  2  -1/2 

%  -  -  u.d  ♦  «0>  .  %  -  -  v-(1  +  V  (3) 
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ua  -  +  U.d  +  o‘) 


-1/2 


2  -1/2 

-  «aUjl  +  o£) 


(M 


in  which  a  and  a  are  the  absolute  values  of  the 
c  a 

slope  at  these  points.  Point  b  is  a  stagnation  point. 

For  aanoothly  curving  foil  the  part  abc  will  be 
represented  in  the  oi-piane  by  a  smooth  continuous  ourve 
passing  through  the  three  known  points  a,  b  and  o.  The 
analysis  will  use  a  circular  arc  passing  through  the  same 
points  as  a  basis  to  perturb  around.  This  choice,  which 
is  admittedly  arbitrary,  leads  to  a  base  flow  for  which 
the  conformal  mappings  are  given  in  terms  of  elementary 
functions . 

3.  METHOD  OF  ANALYSIS.  If  a  relation  between  the 
potential  W  and  the  velocity  uj  can  be  found  by  confor¬ 
mal  mapping,  the  physical  plane  z  can  be  constructed 
from  (2)  by  integration.  The  first  step  of  mapping  is 
the  rotation  of  the  u)-plane  through  an  angle  P  where  P 
is  the  angle  made  by  the  radius  of  the  circular  arc  at 
point  b  and  the  u-axis.  Thus 

<«*  .  u>e-i^  (3) 

* 

It  is  convenient  to  map  the  w  -plane  onto  a  new  half 

*  *  *  * 
plane  w  »  I  +  iT|  .  The  relation  between  the  u»  -  and 
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* 

C  -planes  is  taken  to  be  of  the  form 

*  wc  -  «(C  )  »a 

U)  .  -  i  i  — 

1  -  «(C  ) 


(6) 


where 


In  these  expressions  k  is  the  angle  of  intersection 
between  the  circular  arcs  cla  and  abc. 

The  C  -plane  will  be  as  shown  in  Figure  Id.  The 
circular  arc  abc  in  the  hodograph  plane  will  be 
represented  by  the  line  segment 

0s§‘st,  D  -  0 
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where  the  oonstant  T  is  defined  by 


The  ourve  »bo  which  represents  the  hydrofoil  will 
deviate  from  that  segment  as  shown.  The  free  streamlines 
al  and  ol  will  be  represented  respectively  by  the  por¬ 
tions  s  0  and  T  *  5*  *  *  of  the  {*-axis. 

Nov  define  a  half  plane  C  ■  S  +  iT)  with  the  position 

_  * 

of  points  a,b,e  the  same  as  those  in  the  Q  -plane. 

„  *  „ 

Assume  the  transformation  from  the  Q  -plane  to  the  w- 
plane  to  be  of  the  form 

C*  -  C  ♦  «(C)  (il) 

Here  *(C)  is  an  unknown  complex  analytic  function  of  C 
that  can  be  separated  into  real  and  imaginary  parts  as 
follows 


«(5,T))  -  X({,TU  +  iY(5,T)) 


(12) 


Since  both  C*  and  C  are  real  on  the  lines  al  and 
cl,  then  the  imaginary  part  in  equation  (12)  has  to 
vanish  on  those  lines,  or 

T  S  t  *  • 

Y<*»°>  “  0  ..<(<0  (I3) 
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Thus  Y(5»T1)  is  the  imaginary  part  of  an  analytio 
funotion  and  hence  satisfies  the  differential  equation 


S{ 


+  Yn  "  0 


D  *  0  (l4) 


Assume  that  on  the  part  abo,  Y  will  be  represented  in 
the  general  Fourier  series  form 


Y(?,0)  «  £  [aR  sin(nTrrl)+  bR  oos(nTrr§)] 

n»  0 

0  s  5  s  T  (15) 

The  boundary-value  problem  given  by  equations  (13),  (l4) 
and  (15)  has  the  solution 


r>°° 


Y(5,TU  -  i 


Tr 


)  {a^  sin(nTrrs) 
n-0  n 


b  cos(nTTTs)} 


n 


} 


cos{t( s  -  5  ) } 


-til 

ds  e  1  dt 


The  analyticity  of  c(C)  implies  that 


(16) 


Y 


D 


(17) 
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From  equations  (l6)  and  (17)  the  relation  for  X(£,?l) 
is  found  to  be 


X(S,D) 


1 

n 


s  in  (nT  ns)  +  b 

n 


cos(nTn  a) }  {• 


sin{t( s 


ds  e 


-tT) 


dt 


(18) 


Relations  (l 6)  and  (18)  completely  define  t(C)  in 

terms  of  the  unknown  constants  a  and  b  . 

n  n 

By  the  theorem  of  Schwarz  and  Christoffel,  the 
oomples  potential  is  mapped  onto  the  same  C-plane 
through  the  transformation 

W  -  ce  (C  -  l)2  (19) 

cl 


where  tp  is  the  velocity  potential  at  point  a.  From 
the  definition  of  u>  given  by  (2)  and  (5)  one  finds 
that 


dz  e~lg  1  dW 

3C - 


(20) 
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•  dV 

A  substitution  for  u>  and  for  ^  by  differentiating 
(11)  leads  to  the  integral 


Here 


F(C) 


“c  -  *(C)  \ 

-  r™”i-rn"- 


(22) 


e  (C) 


1  r  (m  -  n) 

*  [Hi  +  ■«"(C)'J 


+ 


(23) 


The  right-hand  side  of  equation  (21)  involves  the 

function  c(C)  which  is  determined  in  terms  of  the 

unknown  constants  a  and  b  .  To  evaluate  the  first 

n  n 

m  of  these  constants  use  can  be  made  of  a  method 
similar  to  that  introduced  by  Naiman  as  summarized  by 
Abbott  and  Von  Doenhoff  (1959)*  Divide  the  arc  abc 
in  both  the  z-  and  C-planes  into  (m  +  1 )  equal 
segments.  Make  a  substitution  for  each  point  at  the 
end  of  every  segment.  For  example,  at  the  end  of  seg¬ 
ment  number  k  the  equation  will  be 
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2|B  -IP  pC  k/ (n+1) 

•*  ^k - V-  J  dt  (24» 

Cc 

k  ■ 

Thus  one  obtains  2(m+l)  equations  involving  the  con¬ 
stants  Sq  f  a^ |  • • • f  t  b^ j  * • • *  b^i  Solving 

these  equations  simultaneously  will  lead  to  the  numeri¬ 
cal  value  of  each  of  these  constants.  The  number  m 
can  be  ohosenas  large  as  the  aoouracy  requires* 

* 

The  determination  of  the  constants  defines  u>  and 
henoe  u>  completely.  The  net  pressure  at  any  point  of 
the  lamina  is  given  by 

P  -  |pu»  {l  -  MS>|2  }  (2 5) 


where  |u)(§)|  is  the  absolute  value  of  u>(§)  computed 
from  (  5). 

The  components  of  the  total  force  on  the  lamina  are 
found  by  integrating  equation  (25 )t 


l»(5)|2  }|f$)« 


(26) 


(27) 
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Thus  the  lift  coefficient  will  be  given  by 


2  ( F  cos  ex  -  F  sin  o) 

„  y  * 

Ct — IS - 


(28) 


Similarly  the  drag  coefficient  can  be  found  from  the 

relation 


2  ( F  cos  a  +  F  sin  a) 

C  .  - - - — g — * - 

P  U_B 


(29) 


4.  HYDROFOILS  WITH  FINITE  CAVITATION  NUMBER.  In  the  case 

that  the  pressure  p  inside  the  cavity  is  different  from 

c 

p^  the  free  streamlines  al  and  cl  seem  to  close  at  a 
finite  distance  and  the  length  of  the  cavity  will  be  no 
longer  infinite.  In  this  case  the  velocity  on  the  free 
streamline  V  will  be  related  to  by  the  relation 


|pv2  -  (ptt  -  Pc)  +  |pU*  (30) 

The  mathematical  model  to  be  used  here  (Figure  2a)  is  the 
same  as  the  one  introduced  by  Wu  (1962)  for  cavity  flow 
around  flat  plates.  Two  curved  plates  dl  and  d'l  are 
supposed  to  extend  to  infinity  where  the  pressure  and 
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velocity  on  both  plates  change  from  p  and  V  at  points 

c 

cl  and  d'  to  ptt  and  as  they  reach  points  1.  In 

general  the  plates  have  a  curvature  that  is  unknown  at  the 
outset.  Hie  complex  potential  (Figure  2b)  and  the  complex 
velocity  planes  are  drawn  with  the  positions  of  points  d 
and  d'  as  indicated.  For  convenience,  however,  in  this 
case  the  complex  velocity  is  defined  as 


U) 


1  d*v  u  .  v 

“  1  dl  "  V  "  1  ? 


(31) 


Define  the  C^-plane  to  be 


c 


1 


i 


* 

Cl) 

a 

-  j 

0)  -  Cl) 

a  c 


*  * 

0)  -  CD 

C 


*  * 

0)  -  U) 

a 


tt/k 


(32) 


where  u>  has  the  same  definition  given  by  (5).  Choose 
the  plates  dl  and  d'l  to  be  of  such  a  shape  to  map  onto 
a  vertical  slit  as  shown  in  Figure  2d.  By  means  of  the 
Schwarz-Christoff el  theorem,  the  C^-plane  can  be  trans¬ 
formed  to  a  half  plane  p  (Figure  2e).  In  general  the 
shape  abc  will  not  transform  into  the  axis  as  indicated, 
but  must  be  mapped  according  to  an  additional  transforma¬ 
tion  similar  to  the  procedure  of  the  foregoing  section. 

The  W-plane  and  p-plane  are  also  easily  related  by  the 
method  of  Schwarz  and  Christoffel. 
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5.  A  RAPID  METHOD  FOR  ESTIMATING  LIFT  AND  DRAG  COEFFICIENTS 


The  method  introduced  here  is  a  means  to  get  a  fast  estimate 
for  the  lift  and  drag  coefficients  based  on  the  assumption 
that  the  free  streamline  for  the  foil  of  interest  is  very 
close  to  that  computed  from  the  basic  flow  for  which  the 
hodograph  boundary  consists  of  two  oiroular  arcs.  In  effect 
a  correction  is  made  on  the  wetted  foil  surface  while  the 
proper  boundary  condition  on  the  free  streamlines  is  ignored. 

First,  assume  that 


U) 


»*  -  e(C) 

c  a 


i  -  g  (i) 


T(C) 


(33) 


on  the  curve  abc  in  the  C~plane  and,  second,  assume  that 


t(S,0) 


00 


(3*) 


From  the  above  two  relations  and  the  definition  of  W  and 
cu  the  physical  plane  is  described  by  the  formula 


2V 


-ip 


U. 


!L ±  »; ( 

J5  F(?)  -WmJg F(6 


{2nmi5/S  } 


) 


dS 


I 


5  (?-U 


^am  exp  {2rrmi§/§c}| 
•00  J 


f3(C) 


d5  +  •  •  • 


By  substituting  for  i  and  z  at  corresponding  points 

the  above  relation  yields  a  set  of  simultaneous  equations 

for  the  complex  constants  a  .  Substitution  of  these  con- 

m 

stants  in  (3*0  will  define  u>  along  abc.  Hie  lift  and 
drag  coefficients  may  then  be  formed  form  (26)  and  (27). 

As  an  example  of  the  method  consider  the  hydrofoil  to 
be  a  circular  arc  subtending  an  angle  n/3  set  in  the 
stream  at  an  angle  a  which  varies  from  0  to  rr/2.  The 
cavitation  number  is  zero.  The  velocity  components  at 
points  a  and  c  will  be 


* 


U) 

a 


+ 


T 


i 


The  values  of  the  constants  X  and  M  are 


\ 


2 
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Table  I  gives  the  values  of  the  constant  N  for  selected 
values  of  ».  The  constants  a^  have  been  evaluated  by 
matching  the  conditions  from  equation  (35)  a,t  ten  points 
along  the  surface  of  the  foil.  The  resulting  values  for 
the  force  coefficients  are  given  in  Table  II.  Figure  3 
shows  that  the  agreement  with  the  computations  of  Rosenhead 
is  very  satisfactory. 
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APPENDIX  A 

PERTURBATION  OF  THE  FLOW  PAST  A  PLAT  PLATE 


In  the  method  of  Tulin  (1953)  the  flow  about  a  super- 
oavitating  foil  is  regarded  as  a  perturbation  of  a  uniform 
flow.  Thus  both  the  attack  angle  and  camber  are  restricted. 

An  alternative  scheme  would  use  the  cavity  flow  about  a  flat 
plate  as  a  basis  for  perturbation.  Such  a  scheme  is  out¬ 
lined  here,  up  to  terms  of  the  second  order.  The  oamber  must 
be  in  some  sense  small,  but  there  is  no  limitation  on  angle 
of  attack. 

A. 1  METHOD  OF  ANALYSIS.  Let  the  cambered  foil  a*b*c*  be 
set  such  that  the  chord  a*c*  makes  an  angle  a  with  the 
direction  of  the  flow  as  shown  in  Figure  4.  Let  a^c*  be 
considered  as  the  z-aiis,  and  the  line  perpendicular  at 
point  a*  be  the  y-axis.  Assume  that  the  maximum  deflection 
of  the  foil  is  h  and  the  chord  length  is  B.  Let  the  small¬ 
ness  parameter  6  be  defined  by 

0  »  h/ B  ■*=  ■<  i 

It  is  now  assumed  that  the  stream  function  ^  satisfies  the 
Laplace  equation 

X,  *%y-°  U.l 
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and  moreover  that  it  oan  be  expressed  in  the  form 

¥(*#y)  -  ♦D(*»y)  +  +  62t2(x,y)  4  ...  (A. 2) 

in  which  tQ(x,y)  is  the  stream  funotion  if  the  foil  is 
replaced  by  a  flat  plate  ao.  Also  let  the  ordinate  y 
for  the  lamina  a*b*c*  be  given  by  the  relation 

y  -  M(x) 

The  first  objeotive  of  the  analysis  is  to  formulate  and  solve 
the  appropriate  boundary  value  problems  for  the  determination 
of  tQ(x,y)  and  the  first-  and  second-order  perturbation 
terms  t^Xjy)  and  t2(x,y)  •  Since  6  is  for  the  moment 
arbitrary,  it  follows  from  (A.l)  and  (A. 2)  that  1Q,  t2» 

...  satisfy  the  Laplace  equation.  Determination  of  the 
boundary  conditions  must  now  be  considered.  Assign  for  the 
streamlines  I*b*a*I*,  I*b*c*I*  the  values  ^  «=  0.  Then  on 
a*b*c*  we  will  have  the  condition 

-  to{*,6Tl(x)]  +  6*1{x,6THx)} 

+  62*2{x,6Ti(x)1  . . 0 

Expanding  the  terms  in  the  right  hand  side  by  means  of  a 
Taylor  series  leads  to  the  relation 

♦.<*•»>  *  ♦«,<*.<»  *  *o„U,0>  ♦  ...  * 

+  6*1(x,0)  +  621Hx)  (*,0) 
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+  •  •  •  + 


0 


♦  fi2*2U,0)  +  ...  - 

By  equating  the  terns  with  like  powers  of  6  we  obtain  the 
relations 

♦o(i,0)  -  0  (A. 3) 

T\(x)  toy(x,0)  4  ^(*,0)  -  0  (A. 4) 

■iyT^U)  iQyy(x,0)  4  T)(x)  tly(x,0)  4  tgU.O)  -  0  (A. 5) 

Let  xQ,yo  represent  the  values  of  x,y  on  the  free 
surface  for  the  case  of  a  flat  plate,  as  indicated  in  Figure  4. 
Assume  that  the  coordinates  x,y  of  the  free  streamline  for 
the  cambered  foil  a*b*c*  to  be  given  in  the  form 

y  -  yG, 

x  m  xo  *  +  &252(y)  +  •••  (a. 6) 

Using  the  fact  that  the  free  surface  is  a  streamline,  we 
obtain  the  relation 
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¥t«.  +  6?1(y)  +  62?2(y)  +  ...,yol  -  t0{xo  +  ti1(y)  +  62§2(y)  +  ...,yo) 
+  6*ltxo  +  4Sl^y  ^  +  62?2^  4  *  *  *  »yo^  +  ®2*2^Xo  +  6^1(y )  ♦  62$2(y)  ♦  . . .  ,yl 


+  ...  -  0 


An  expansion  by  Taylor  series  as  before  leads  to  the  following 
boundary  conditions 

*oUo'yo)  ’  0  (A*7) 

§l(y)  to*(Vyo)  +  ♦l(xo’yo)  -  0  (A*8) 

?2(y)  *ox(xo’yo)  +¥l{ y)  ♦oxx(xo»yo) 

+  5l(y)*iy(*o»yo)  +  ♦2(xo»yo)  ■  0  (A‘9) 

One  more  condition  on  the  free  surface  is  that  the  resultant 

velocity  V  is  constant,  or 
c 
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+  *2?2(y)  +  .. 


dx 


tx  +  6§1(y)  +  62§2(y)  +  ...,y} 


3y 


A  substitution  for 


and 


a¥/ay 


from  (A. 2)  leads  to 


the  following  results « 


♦o,x<V>'„>*6{l<>'>  *,>x.<V’'o)**%><’'> 


*  I*25!*’’1  *xxx<VV*  ••■*  (*lx(V>,o>*  ‘2ht,’>  *1X1  •’W 


+  62* 


2x 


2 


*oyUo’yo)  +  *Vy)  *oyx(xo’yo)  +  6  52(y)  ♦oy*Uo’yo) 


|«251(y) 


*oyxx(xo’yo)  +  "‘  +  6*ly(Vyo)  +  *2?i(y) 


^lyx (xo*y0) 


+  fl^2y(xo'yo)  + 


(A. 10) 
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Equating  terms  with  like  powers  of  6  as  before  results  in 
the  relations 


Vc  ‘  (*o*)2  +  <*oy>2  ‘A*11 

Vy)(*ox*oxx  *  ♦0y*oyx)  +  *ox*lx  +  *oy*ly  “  0  (A.12) 

§l(y)*oxx  +  *lx  +  2?2<y)*ox*oxx  +  ?l(y)*ox*oxxx  + 

*  25l(y)*ox*lxx  +  2*ox*2x+5l<y>*oyx  + 

+  252(y)toytoyx  +  (y  ^oy^oyxx  + 

+  25l(y)*oy*lyx  +  2toy*2y  "  0  (A*1?) 


In  these  expressions  the  derivatives  of  the  stream  function 
terms  such  as  tox»  t1:t>  •••»  are  to  be  evaluated  on  the 
unperturbed  free  streamline  J0»yo*  Conditions  (A. 3), 

(A. 7)  and  (A. 13)  together  with  the  differential  equation 


*0X31  +  *oyy 


(A. 14) 


formulate  the  boundary  value  problem  for  the  zero  order  term 
♦o(x,y).  The  solution  for  this  problem  for  zero  cavitation 
number  was  given  by  Rayleigh  (see  Lamb,  1932)  using  the 
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hodograph  method.  Since  we  are  going  to  use  this  solution 
in  a  different  form  for  the  development  of  the  first-  and 
second-order  terms,  we  will  present  it  in  the  following 
section. 

A.  2  SOLUTION  TOR  THE  2ER0  ORDER  TERM.  The  physical  plane 
z,  the  complex  potential  plane  W  and  the  dimensionless 
complex  velocity  plane  u>  have  the  same  definition  as 
before  and  are  shown  in  Figure  5  . 

The  new  plane  W-*  is  also  shown  in  the  same  figure. 
Define  the  half  plane  t  to  be 


t  -  -(l/u)  +  u>)  (A. 15) 

By  using  the  Schwarz- Christoff el  theorem  the  W  1  is  mapped 
onto  the  half  plane  t  and  thus  the  relation  between  t  and 
W  is  given  by  the  formula 

C 

w  =  - *  ( A.  l6 ) 

(t  -  2  cos  o / 


where  C  is  a  constant  dependent  on  the  dimensions  of  the 
lamina.  A  substitution  in  (A. 15)  with  the  value  of  to  and 
use  of  Eq.  (A.l6)  leads  to  the  relation 


t 


-  2  C 


U  (t  -  2  cos  cv) 


dt 

"5  dz 


ujt  - 


q 

2  cos  a) 
2C 


dz 

dt 


or  we  get  the  following  differential  equation  to  relate  the 
z  and  t  planes t 
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7 


d  .2  2Ct  /drl 

dt'  U  (t  -  2  cos  a dt  U^(t  -  2  cos  Of)^ 


0  (A. 17) 


The  above  has  the  following  solution i 


z 


_C_ 

U_ 


P  t 


°2 


-t  +  V t2 -  4 
(t-  2  cos  «)■ 


dt 


where  the  complex  constant  C  is  determined  by 


(A. 18) 


1 

C 


-1 


BU 


00 


p2 


-  2 


-t  +  -  4 

(t-  2  cos  «)^ 


(A. 19) 


Substitution  for  t  from  (A15)  results  in  the  integral 


z 


r>  O) 


_C_ 

U„ 


6  4 

w  -  w 


2 

u>  +  1 


|_j  u>{ -  (2  cos  o')  0)  +  l}- 


du> 


(A. 20) 


Similarly  by  substitution  for  t  in  (A.l6)  we  find  W  in 
terms  of  u),  thus 


W 


Cuj 

{«  +  (2  cos  Of)  <o  +  1} 


(A. 21) 
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Since  W  -  <pQ  +  itQ  ,  one  sees  that 


U) 


♦o(x,y)  -  Img  — £ 


u>  +  (2  cos  «)  a)  +  1, 


( A .  2  2 ) 


The  relations  (A. 20)  and  (A. 21)  give  t  and  V  in  terms 
of  the  parameter  <u  and  hence  a  complete  solution  for  the 
cero  order  term  is  accomplished. 


A. 3  SOLUTION  FOR  HIGHER  ORDER  TERMS.  Consider  now  the 
first  order  term  which  is  represented  by  the  following 

boundary-value  problem i 


^lxx  +  * lyy 


-  0 


(A. 23a) 


T){x)  tly(x,0)  +  ^(x.O)  -  0 


(A. 23b) 


Vy)  *ox(xc’yo)  +  *l(xo’yo>  "  0 

5l(y,[toxUo*yo)  *on(lo’3'o)  +  ♦oy(xo'yo)  *oyx (i  o  ’yo  >]  - 
♦  ♦ox(xo»yo)  *lS(xo’yo>  +  ♦oy<*o*o)  ’  ly ^xe  'o)  "  0 


{A. 23c) 


(  .23d) 
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The  last  two  conditions  can  be  combined  into  the  condition i 


[♦ox(xo’yo)  ♦o**uo’yo)  ♦  *y 
♦♦ox(xo»3ro,Kx(xo**o)  *lx(xo'yo)  +  ♦oy 


(xo’yo>  ♦ox*<xo»yo)] 
(xo*yo)  ♦ly(xo'yo)]" 


0 

(A. 24) 


Due  to  the  complexity  of  the  s-plane,  we  will  use  Eq.  (A. 20) 
to  map  our  domain  into  the  semi-circular  hodograph  planet 


U>  -  +  lu>2  .  U  +  iV 

The  last  notation  differs  from  that  of  the  main  section  of 
the  report.  In  irfiat  follows  it  is  convenient  to  introduce 
the  abbreviations  listed  belowt 


P  - 


t  u 
ou  x 


+  t  v 

ov  X 


*ouuy 


+  ♦  V 
ov  y 


-  *luux 


+  *lvTx 


-  *luuy 


+  ♦lvTy 
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The  previous  boundary  value  problem  now  assumes  the  formi 


*luu  +  ^lvv 


(A. 25a) 


Tl(u,0) 


♦lu(u'0)  I?  +  *lv(u’0)  f? 


+  -  0 


(A. 25b) 


?1(u,v)P  ♦  t1(u,v)  -  0 


(A. 250) 


51(u,v)<;  p 


Puux  +  Vx 


+  Q  Q  u  +  Q  v 
I  u  y  v  y 


+  PR  +  QR  >  ■  0  (A.25d) 


The  last  two  boundary  conditions  are  valid  on  the  circum- 
ference  of  the  semi-circle  in  the  hodograph  plane  and,  as 
before,  can  be  combined  into  a  single  formulas 


r  , 

(u,v)  j  P 
A  I 


P  u  +  P  v 

u  X  V  X 


r 

QiPuuy  +  Vy 


l  +  P  [PR  +  QR]  -  0 


(A. 26) 


The  derivatives  u  ,  v  ,  u  ,  v  are  determined  from  Eq. 

*  "  7  7 

(A. 20).  It  is  obvious  that  this  linear  boundary  value 
problem  can  be  solved  numerically  in  the  hodograph  plane  by 
taking  as  small  a  mesh  as  the  accuracy  requires. 
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By  a  similar  procedure  the  formulas  for  the  second- 


order  term  tg(x,y)  are  developed  and  are  found  to  bet 


*2uu  +  *2vv  "  0 


(A, 


■|ll2(u,v)  [uyQu  +  +  T)(u, v)S  +  ♦  (u,0)  -  0 

52(u,v)P  +  ■|5^(u,v)[u3[Pu  +  vxpv]  +  (Uf v)Q  +  *2(u,v)  -  0 

5l(u,T,[“lPu  ♦  Vv]  *  r2  ♦  2'2<U'V)PKPU  *  %p.]  * 

*  Sl<“.'>pK<“IPu  *  Wu  ♦  *  Vv,vJ  * 

*2S1(u,»)p[uJPu  *  »,PV]  •  2P[V2U  *  V2.l  ♦ 

2 

+  §l(u,v)[uyPu  +  vyPv]  ♦  2§2(u,v)Q[uyPu  +  vyPv]  + 

+  5l(u’v)<5[uy(u*Pu  +  vyPv)u  +  vy(uxPu  +  vxPv)v]  + 
+2C1(u,v)pJuyRu  +  vyRy]  +  2Q[uyt2u  4  vyt2y]  -  0 


(A 


(A, 


(A, 


In  principle  the  above  problem  can  also  be  done  numerically 
in  the  semi-circular  hodograph  plane.  The  equations  for  the 
interior  points  of  the  mesh  will  be  the  same  but  those  for 
the  boundary  points  will  have  to  satisfy  the  conditions  in 
(A. 27)  instead  of  (A. 23). 


27a) 

.27b) 

2?  c) 


27d) 


28 


The  relation  for  the  free  streamline  is  given  in  (A. 6) 
where  5^(y)  and  5g(y)  are  to  be  determined  from  Eqs. 

(A. 8)  and  (A. 9).  The  actual  velocity  components  U  and  V 
at  any  point  are  computed  from  the  formulas  t 


/»♦  0 
-N* 


+ 


The  local  pressure  at  any  point  on  the  lamina  may  then  be 
computed  by  use  of  Bernoulli's  equation,  and  the  total  foroe 
obtained  by  an  integration. 
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APPENDIX  B 
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LIST  OF  SYMBOLS 

Fourier  coefficients 

Chord  length 
Drag  coefficient 

Lift  coefficient 

Function  defined  by  (22) 

Force  components 

Function  defined  by  (23) 

Maximum  camber  of  foil 
Constant  defined  by  (9) 

Constant  defined  by  (10) 

Pressure  inside  the  cavity 

Pressure  at  infinity 

Complex  variable 

Constant  defined  by  T  =»  1  -  M/N 
Velocity  magnitude  at  infinity 

Velocity  components 

Velocity  magnitude  on  the  free  streamline 
Complex  potential  W  =  cp  +  if 
Cartesian  coordinates 
Physical  plane,  z  •  x  +  iy 
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Angle  of  attack 

Angle  through  which  hodograph  ia  rotated 
Smallness  parameter,  6  -  H/B 
Complex  variable  e ( C )  -  X  +  iY 
Complex  variable,  Q  *  §  +  ill 
Complex  variable 

Angle  of  intersection  of  circle  arcs  in  the 
approximate  hodograph  plane 

Constant  defined  by  (8) 

Complex  variable 

Fluid  density 

Cavitation  number 

Function  defined  by  (3*0 

Velocity  potential 

Stream  function 

Dimensionless  complex  velocity  defined  by  (2) 
Complex  variable  defined  by  (5) 
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TABLE  I. 

VALUES  OF  N 

Of 

N 

0° 

0.19253 

o 

o 

0.49074 

60° 

1.54024 

90° 

IO.OO725 

TABLE  II.  LIFT  AND  DRAG  COEFFICIENTS 


Attack 

Angle 

« 

Lift  Coefficient 

Drag  Coefficient  C^ 

Present 

Theory 

Rosenhead 

Present 

Theory 

Rosenhead 

0 

0.8081 

— 

0.0406 

— 

O 

O 

0^ 

.7091 

0.7042 

.4109 

0.4038 

O 

O 

— 

.6620 

— 

.5484 

O 

O 

— 

.5784 

— 

.6796 

6o° 

.5132 

— 

.8620 

— 

70° 

— 

.3226 

.8740 

90° 

0 

0 

- 

.9531 

_ 

.9438 

_ 
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TABLE  III.  VALUES  OF  THE  COEFFICIENTS  a 

_ in 


m 

o 

o 

■ 

a-30° 

a.6o° 

a.  90° 

-4 

.07707 

.03602 

-.00606 

-.00102 

-3 

.02856 

-.19926 

-.05532 

-.03836 

-2 

.13751 

-.03389 

.00209 

.00225 

-1 

.17162 

-.16834 

-.03658 

-.05085 

0 

-.23252 

-.07860 

.00392 

.00428 

1 

.054419 

-.20105 

-.05209 

-.03745 

2 

-.030412 

-.08845 

.00541 

.00302 

3 

-.116138 

-.22033 

-.04130 

-.05387 

4 

.04903 

-.07925 

.00384 

.00117 

5 

-.07337 

-.11827 

-.05260 

-.04211 
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1.  Conformal  mapping  planes  for  cambered  hydrofoil 
at  tero  cavitation  number. 
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Figure  2.  Conformal  mapping  planes  for  hydrofoil  at  finite 


cavitation  number. 
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